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Abstract: 

We have computed the cross section for tt production by gluon-gluon fusion at next-to- 
leading order (NLO) in the QCD coupling, keeping the full dependence on the tt spins. 
Furthermore we have determined to the same order the spin dependent cross sections for 
the processes g + q(q) —*ti+q(q). Together with our previous results [JT]] for q + q ^ tt(g) 
these results allow for predictions, at NLO QCD, of the hadronic production of tt pairs 
in a general spin configuration. As an application we have determined the degree of 
correlation of the t and t spins at NLO, using various spin quantisation axes. 
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In this Letter we report on the calculation of the cross section for tt production by gluon- 
gluon and by (anti)quark-gluon fusion at order a, in the QCD coupling, keeping the full 
information on the spin state of the tt system. The results presented in this work constitute 
the last missing ingredient to analyse, at NLO QCD, tt production including spin effects 
at hadron colliders. Together with the corresponding results for tt production by quark- 
antiquark annihilation |JlJ], they extend previous results for the spin-averaged differential 
tt cross section [@, g, g, |5p . 

Analysis of top quark spin phenomena at the upgraded Tevatron and at the Large 
Hadron Collider (LHC) will play an important role in investigating the interactions of 
these quarks. (For a recent overview of the perspectives of top quark spin physics at 
hadron colliders and an extensive list of references, see, e.g., ref. In the following we 
first describe the salient features of our computation. Then we determine the correlation 
of the t and t spins at NLO for a number of spin bases. 

The theoretical description of tt production by proton-proton and proton- antipro ton 
collisions at next-to-leading order in the QCD coupling amounts to considering the fol- 
lowing parton reactions: quark- antiquark annihilation 

q(pi) +q(p2) -> t(kt) +t(k), (i) 
q{pi) +q(p2) -> t(k t ) +t(k T ) + g (k 3 ), (2) 

which is the dominant production mechanism at the Tevatron, gluon-gluon fusion 

g(pi)+g(P2)^t(k t )+t(ki), (3) 

g(pi) +g(Pi) -> t(k t ) +t{k) +g(k 3 ), (4) 
which dominates at the LHC, and 

g(pi) +q{pi) -> t(k t ) +f(kt) +q(h), (5) 

g(Pl) +q(Pl) -> t(k t ) +t(k) +q(k 3 ). (6) 

In the Standard Model the main top decay modes are t — > bW — > bqq'.biVi. Among 
these final states the charged leptons, or the jets from quarks of weak isospin —1/2 orig- 
inating from the W decay, are the most powerful analysers of the polarisation of the top 
quark. A complete next-to-leading order QCD analysis thus consists of treating the parton 
reactions 

gg,qq^tt^bb + 4f, (7) 

gg,qq^tt^bb + Af + g, (8) 

g + q(q)^ti^bb + 4f + q(q), (9) 

where f = q, i,V£. To leading order in a s only the reactions (|7|) contribute. In view of the 
fact that the total width T t of the top quark is much smaller than its mass, T t /m t = 0(1%), 
one may analyse the above reactions using the so-called narrow width or leading pole ap- 
proximation [g, |J]. In the case at hand this approach consists, for a given process, in an 
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expansion of the amplitude around the poles of the unstable t and t quarks, which corre- 
sponds to an expansion in powers of T t /m t . Only the leading term of this expansion, i.e., 
the residue of the double poles is considered here. In this framework the radiative correc- 
tions to (0) can be classified into so-called factorisable and non-factorisable corrections. 
In the case of the gluon radiation processes (§) already the lowest order contributions 
to the respective squared matrix element 1 5Vf | 2 decompose into these two groups. We 
compute the factorisable corrections^ for which the squared matrix element M is of the 
form 

\M | 2 oc Tr [pflp] = pa'a^aa',pp'Pp'p- (10) 

Here R denotes the density matrix for the production of on-shell tt pairs and p,p are 
the density matrices describing the decay of polarised t and t quarks, respectively, into 
specific final states. The subscripts in eq. ([TT]) denote the t, t spin indices. Because the 
colours of the final state partons are summed over, p,p are unit matrices in the colour 
indices of t and t. This implies that the colour indices of t and t are summed over in the 
computation of R. Note that both the production and decay density matrices are gauge 
invariant. The one-loop QCD corrections to the density matrices of semileptonic decays 
of polarised top quarks and of t — > W + b can be obtained from the results of refs. [ |T0[ ] 
and [ p] , [T2"| ], respectively. 

In the calculation of the production density matrices for the reactions (Q) - (0) we 
follow our approach [|TJ to computing those of reaction ([TJ) and reaction (Q). The density 
matrix for the production process @ is defined in terms of the transition matrix element 
as follows: 

R aa>,W = W- E (taW\gg}(gg\^\ta>tv), (11) 

ly 88 colors 

initial spins 

where the factor N gg = 256 averages over the spins and colours of the initial pair of gluons. 
The matrix structure of R^ is 

+ Clf(&)aa>(0 j )w, (12) 

where o' are the Pauli matrices. Using rotational invariance the 'structure functions' 
B (gg)^(gg) and c te) can be f urt her decomposed. The function A^s) = Tr [R^}/4, de- 
termines the gg — > tt differential cross section with tt spins summed over. Because of 
parity invariance the vectors B^ 88 \B^ can have, within QCD, only a component nor- 
mal to the scattering plane. This component, which amounts to a normal polarisation of 
the t and t quarks, is induced by the absorptive part of the scattering amplitude which is 
non-zero to order a^. The normal polarisation is quite small, both for tt production at the 

Tevatron and at the LHC Jl3| , [14| ] . The functions c\ 88 ^ encode the correlation between the 
t and t spins. The production density matrices R\88>g) , R(g4>Q) , and R(g9>9) for the reactions 
3 The non-factorisable NLO QCD contributions were studied for gg and qq initial states in ref. [p|] . 
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(Q), (@), and @, can be defined and decomposed in an analogous fashion, with the colour 
and spin degrees of freedom of the final state gluon, quark q, or antiquark q being summed 
over. 

In Born approximation R^ gg > was given, e.g., in ref. [JTSj]. In the computation of the 



order oc^ contributions we used dimensional regularization to treat both the ultraviolet 
and the infrared/collinear singularities. Renormalisation was performed using the MS 
prescription for the QCD coupling a s and the on- shell definition of the top mass m t . The 
remaining soft and collinear singularities in R( gg ) , which appear as single and double poles 
in e = (4 — D)/2, are cancelled after including the contributions of R( gg > g ) in the soft and 
collinear limits and after mass factorisation. For the latter we used the MS factorisation 
scheme. The soft and collinear singularities of R^88,g) were extracted by employing a 
simplified version of the phase-space slicing technique (cf. ref. [[Ij] for details). 

As a check of our calculation we computed the total cross section for gg — > tt+X at 
NLO. If one identifies the MS renormalisation scale /j with the mass factorisation scale /jp 
and neglects all quark masses except for m t , then one can express the parton cross section 
in terms of dimensionless scaling functions [0]: 

d gg (s,mj) = 4[4 0) 0l) +4;ra,(/ij ) (ri) +ff g ) (T ] )ln( P 2 /mt))], (13) 
m t 

where s is the parton cm. energy squared and 

" = ^- 1 - <14) 

Our results for these functions are shown in Fig. 1. The dotted line is the Born result 

fgg (t|), the full line shows the function fgj, (tj), and the dashed line is fgg We 

find excellent agreement with the corresponding results for fgg^\v\) of ref. [Q, |]] and, 
after accounting for the different renormalisation schemes used by us and in ref. 0], with 

The determination of R^si^) to order oc^, describing the (anti)quark-gluon subpro- 
cesses (||) and (§), involves initial state collinear singularities that we mass-factorised 
using the MS factorisation scheme. The parton cross section d gq (s,mj) takes the form 



$ Sq = +/ff(Tl)ln(A/ 2 M 2 ))] • (15) 

mf 

The functions fgq (tj) (full line) and fgq (tj) (dashed line) are shown in Fig. 2. They agree 
perfectly with the corresponding results of ref. [0, |5p . The density matrix can be 

obtained from Rww by exploiting the charge-conjugation invariance of QCD. 

In order to exhibit the tt spin-correlation effects we consider now the following set of 
observables: 

Oi=4s r s f , (16) 
2 = 4(k f -s,)(k F -s F ), (17) 
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Figure 1: Dimensionless scaling functions fgg\r\) (dotted), fgl\r\) (full), and fgl\v\) 
(dashed) that determine parton cross section d gg . 
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Figure 2: Dimensionless scaling functions fgq(r\) (full) and fgq(r\) (dashed) that deter- 



mine a 
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3 = 4(pi.s,)(pi-s f -), (18) 

04 = 4(pI. Sj )( P f- Sf ), (19) 

where s f = (a<g> 11) /2 and sj = (11® g)/2 are the t and t spin operators and the factor 
4 is conventional. The expectation values of the observables 0%, O3 and O4 determine 
the tt spin correlations using different spin quantisation axes. We consider here the fully 
inclusive case, i.e. integrate over the appropriate full final state phase space. Observable 
02 corresponds to a correlation of the t and t spins in the helicity basis, while O3 correlates 



the spins projected along the beam line in the parton c.m.s. The 'beam-line basis' [ |16| ] 
used in O4 refers to spin axes being parallel to P2 and pj*, respectively. Here P2 (pj*) 
is the direction of the downward (upward) collider beam in the t (f) rest frame which we 
define by performing a rotation-free boost from the parton c. m. frame. The expectation 
values of 02,03, and O4 are related to the perhaps more familiar double spin asymmetries 
of the cross section as follows: 

4((a SfKb S ' )} - o(TT)+o(U)+o(U)+o(iT)' ( ) 

where the arrows refer to the spin state of the top quark (antiquark) using the unit vector 
a (b) as quantisation axis. 

For the gg initial state we write the unnormalised expectation value of a spin-correlation 
observable in analogy to eq. (P~3|) as follows: 



6 gg (o) gg = 4[s£ } 0l) +47m s (g«(Ti) +^ ) (Ti)ln( A , 2 M 2 ))] , (21) 
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Our results for the functions g g °J (r\) , g^g (r\) and ggg (x\) are shown for the four observ- 
ables O1-O4 in Figs. 3-6. In each figure, the dotted line is the Born result g g °J (r\), the full 

line shows the function g g l J (r|), and the dashed line is gg]) (r|). 

The QCD corrections to the unnormalised spin correlations are large close to thresh- 
old. This behaviour is due to the factor 6 gg , which, in this order of perturbation theory, is 
non-zero at threshold due to Coulomb attraction. The corrections to the normalised par- 
tonic expectation values of all spin observables are tiny for small T|. This is also the case 
for the contributions from qq annihilation [[]]]. For larger T], the spin correlations show a 
rich structure and the QCD corrections to the normalised expectation values can become 
significant. For hadronic observables, these effects are damped by the parton distribution 



functions and the QCD corrections are of the order of a few percent QT7|]. 

For the quark-gluon process @ we represent the unnormalised expectation value of 
an observable in the form 

a gq (0) gq = ^fcgft) +|^ ) (ri)ln( A / 2 M 2 ))] . (22) 
mf 

We have plotted in Figs. 7,8 the functions ggq\r\), ggq\^) corresponding to the two 
observables O2, O3. The contributions of this channel to the spin correlations are very 
small. 
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Figure 7: Dimensionless scaling functions ggq(T\) (full) and ggq(r\) (dashed) that deter- 
mine the expectation value o g q(02) gq . 




In summary, we have computed the spin density matrices describing tt production 
by gluon-gluon and by (anti)quark-gluon fusion to order a^. Together with our previous 
result on R qq at NLO this work now provides the ingredients for a complete description of 
spin effects in the hadronic production of top quark pairs at NLO in the strong coupling, 
allowing for a more precise investigation of the tt production and decay processes. 
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